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(57) ABSTRACT 

The stability of a recursive estimator process (e.g., a Kalman 
filter is assured for long time periods by periodically reset- 
ting an error covariance P(t n ) of the system to a predeter- 
mined reset value P r . The recursive process is thus repeti- 
tively forced to start from a selected covariance and continue 
for a time period that is short compared to the system’s total 
operational time period. The time period in which the 
process must maintain its numerical stability is significantly 
reduced as is the demand on the system’s numerical stability. 
The process stability for an extended operational time period 
T 0 is verified by performing the resetting step at the end of 
at least one reset time period T r whose duration is less than 
the operational time period T 0 and then confirming stability 
of the process over the reset time period T r . Because the 
recursive process starts from a selected covariance at the 
beginning of each reset time period T r , confirming stability 
of the process over at least one reset time period substan- 
tially confirms stability over the longer operational time 
period T 0 . 

30 Claims, 6 Drawing Sheets 
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STABLE AND VERIFIABLE STATE 
ESTIMATION METHODS AND SYSTEMS 
WITH SPACECRAFT APPLICATIONS 

GOVERNMENT RIGHTS 5 

The invention described herein was made in the perfo- 
rmance of work under NASA contract number NAS5-98069 
and is subject to the provisions of Section 305 of the 
National Aeronautics and Space Act of 1958 (72 Stat. 435; 

42 USC 2457). 10 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates generally to state estimation 15 
processes and, more particularly, to spacecraft subsystems 
which use such processes. 

2. Description of the Related Art 

Kalman filtering is an estimation technique that combines 
a knowledge of the statistical nature of a system’s measure- 20 
ment errors with a knowledge of the system’s dynamics (as 
represented, for example, in a state space model) to arrive at 
an estimate of the system’s state. In particular, a Kalman 
filter combines a current measurement y(t n ) of a system state 
parameter x(t n ) with measurement and state predictions 25 
y*(t„-) and x*(t n “) of the parameter x(t n ) that are based on 
past measurements to thereby provide a filtered estimate 
x*(0 of the parameter x(t n ). As indicated by the time term 
t„, the filter successively and recursively combines the 
measurements and predictions to obtain estimates with a 
reduced variance (wherein t~ and t n + refer respectively to 
times just before and after each time t n ). 

In an exemplary spacecraft application, a Kalman filter 
processes attitude and inertial rate measurements to provide 35 
an attitude estimate. Although the filter is known to be an 
especially effective state estimator for these systems, the 
process tends to have numerical instabilities when it is 
implemented with digital processors and continued over 
extensive time periods. The instability has a source in the 4Q 
finite accuracy of digital processing and it represents a 
danger to long-term spacecraft missions that repeat the same 
operational mode for years (e.g., as in communication 
processes of geosynchronous spacecraft). In response to this 
danger, computationally expensive filter formulations (e.g., 45 
Joseph formulation and UDU factorization) have been used 
to increase numerical stability. These approaches, however, 
increase the use of processing time which is a limited and 
tightly budgeted spacecraft resource. 

In addition, these stability solutions have a further prob- 50 
lem. Because of the high cost of spacecraft missions and the 
critical nature of spacecraft attitude control systems, it is 
especially desirable to verify the stability of all such pro- 
cesses over time periods that correspond to those in which 
the mission will be operational. When these operational time 55 
periods are measured in years, such testing is not feasible 
and, thus, these solutions are not verifiable. 

SUMMARY OF THE INVENTION 

The present invention is directed to computationally 60 
simple methods for enhancing and verifying the long-term 
stability of a recursive process that generates estimates of a 
system’s state. The stability is enhanced by periodically 
resetting an error covariance P(t n ) of the system to a reset 
value P r that exceeds the operational values P 0 of the 65 
process. The recursive process is thus repetitively forced to 
start from a selected covariance and continue for a time 


2 

period that is short compared to the system’s total opera- 
tional time period. Because the time period over which the 
process must maintain its numerical stability is significantly 
reduced, the demand on its inherent stability is also signifi- 
cantly reduced. 

The stability for an extended operational time period T a is 
verified by performing the resetting step at the end of at least 
one reset time period T r whose duration is less than the 
operational time period T 0 and then confirming stability of 
the process over the reset time period T r . Because the 
process starts from a selected covariance at the beginning of 
each reset time period T r , confirming stability of the process 
over at least one reset time period substantially confirms 
stability over the longer operational time period T 0 . 

Spacecraft system embodiments of the invention are also 
provided. 

The novel features of the invention are set forth with 
particularity in the appended claims. The invention will be 
best understood from the following description when read in 
conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIGS. 1A and IB are flow charts that illustrate method 
embodiments of the present invention for respectively 
enhancing and verifying the long-term stability of a recur- 
sive estimator process; 

FIG. 2 is block diagram that shows an exemplary recur- 
sive estimator process; 

FIG. 3 is a diagram that illustrates the methods of FIGS. 
1A and IB as they relate to an error covariance of the 
recursive estimator process; 

FIG. 4 is a block diagram of a spacecraft attitude control 
system for practicing the methods of FIGS. 1A and IB; and 

FIG. 5 is a diagram of a spacecraft that includes the 
attitude control system of FIG. 4. 

DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENTS 

The present invention is directed to methods and struc- 
tures that enhance and verify the long-term stability of 
recursive estimator processes, e,g., Kalman filter processes. 
In particular, the flow chart 20 of FIG. 1A recites process 
steps 22-26 of a method embodiment that enhances the 
long-term stability and the flow chart 30 of FIG. IB recites 
process steps 22-25, and 32-34 of a method embodiment 
that verifies the long-term stability. Because these flow 
charts contain references to recursive process terminology, 
however, it will facilitate the description and understanding 
of their methods to precede that description with the fol- 
lowing review of Kalman filter processes. 

FIG. 2 illustrates a process 40 that is based on an estimate 
update equation 

x*(C)=x*(t n -) + k(Q{ y(0-y*(V-)}, (1) 

in which a state prediction x*(t„“) is updated by a portion 
k(t n ) of a residue which is the difference between a mea- 
surement y(t„) and a measurement prediction y*(t„“) to 
thereby form an estimate x*(t„ + ). The portion k(t„) is typi- 
cally referred to as the Kalman gain and is calculated as 
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P(Q=E{(err X){err X) T }. 


( 6 ) 


k{t n ) = 


( 2 ) 


crf(r-) + cr2 


o?(tn) = C 1 -k(t n ))o^*(t n ). 


(3) 20 


These processes are expressed in a more general matrix 
form in FIG. 2 where they are partitioned into an estimator 
42 and a gain calculator 44 that supplies a gain K(t„) to a 25 
multiplier 46 of the estimator. The estimator 42 receives a 
measurement Y(t n ) at an input port 48 and provides an 
estimate X*(t n + ) for a time just after the measurement to an 
output port 50. From this estimate matrix, a state prediction 
X*(t„“) is formed for a time just before the next measure- 30 
ment and this state prediction is provided to a summer 52. 

Because the state of the system typically varies dynami- 
cally between measurements, the estimate X*(t n + ) must be 
extrapolated over time to form the state prediction X*(t n _). 
This extrapolation is accomplished by passing the estimate 35 
X*(t„ + ) through a delay 54 and multiplying it by a state 
transition 56 which contains extrapolation information in the 
form of a state transition O^t^). The state prediction 
X*(t„“) is also multiplied by a measurement converter 58 to 
form a measurement prediction Y*(t n “) which is provided to 40 
a differencer 60 where it is differenced with the measure- 
ment Y(t n ). The measurement converter H(t n ) conditions the 
state prediction X*(t n “) so that its elements correspond to 
those of the measurement Y(t n ) and can be properly differ- 
enced with it. 45 

The dififerencer 60, therefore, generates a residue Y(t n )- 
H(t„)Y*(t n “) which is then multiplied in the multiplier 46 to 
form a correction K(t„){Y(t n )-H(t n )Y(t n - )} that will be used 
to update the estimate. The updating is performed in the 
summer 52 where the correction is summed with the state 50 
prediction X*(t n “) to generate the updated estimate 


X*(t /J + )=X*(r„-)+^(t„){Y(t„)-H(t„)Y*(t„-) 


(4) 


p(r„)=£{(Y(o-Y*(t„))(Y(g-Y*(g) r } 


(5) 


When these replacements are made in equation (2), the 
calculation of gain becomes 


in which cr** 2 (t n _ ) is the estimate variance (i.e., uncertainty 
of the estimate) and a m 2 is the measurement variance (i.e., 
uncertainty of the measurement). 

As the process is recursively continued, the estimate 
variance is reduced below the measurement variance so that 10 
the gain k(t n ) declines to a value much less than one. It is 
apparent from equation (1), therefore, that a large portion of 
the residue {y(t n )-y*(t n - )} is initially used to update the 
state prediction x*(t„“) into the updated estimate x*(t n + ) but 
this portion decreases as the process continues (i.e., the 15 
weighting given to new measurements is successively 
reduced). In particular, the estimate variance is reduced at 
each measurement y(t n ) in accordance with 


y+i?(or 


(7) 


in which the the error covariance P(t n ) is appropriately 
modified by the measurement converter H(t„) and its trans- 
form H r (t„). The extrapolation process of equation (7) is 
performed by the gain calculation block 62. 

The error covariance P(t„) is updated between measure- 
ments in an update block 64 wherein the estimate variance 
a/ 2 (t n “) of equation (3) has been replaced by the error 
covariance P(t„) as modified by the measurement converter 
H(t„). Accordingly, the updated error covariance P(t„) is 
processed as 




( 8 ) 


The updated error covariance is then passed through delay 
65 and extrapolated over time to form an error covariance 
P(t„“). As shown in extrapolation block 66, this extrapola- 
tion is processed as 




(9) 


at the output port 50. 

The gain calculator 44 performs similar updating and 55 
extrapolation processes for the gain K(t n ). The estimate 
variance cr** 2 (t n _ ) and the measurement variance a m 2 of 
equation (2) are respectively replaced by an error covariance 
P(t w ) and a measurement-noise covariance R(t n ) which are 
measures of estimate error and of measurement error. In 60 
particular, the error covariance P(t n ) is defined as 


in which E is expectation and T is the transpose operation. 

If the difference between a state parameter and its estimate 65 
is referred to as estimation error and shortened to err, then 
equation (5) becomes 


and is effected with the state transition O^t^) and its 
transform. In addition, plant or process noise, which 
includes modeling errors as well as actual noise and system 
disturbances, is summarized in a process-noise covariance 

In summary, the gain calculator 44 updates the error 
covariance P(t„) in update block 64, delays it in delay 65 and 
extrapolates it in extrapolation block 66 to generate the error 
covariance P(t„“). The updated, delayed and extrapolated 
error covariance in extrapolation block 66 is then combined 
with the measurement-noise covariance R(t„) in the gain 
calculation block 62 to calculate a gain K(t„) which is 
provided to the estimator 22 for processing the current 
measurement Y(t„) at the input port 28 into an updated 
estimate X*(t M + ) at the output port 50. The update block 64 
indicates that the error covariance P(t n ) will reduce over 
time to steady state operational values P 0 . 

Having reviewed the background recursive filtering 
terminology, attention is now redirected to the flow chart 20 
of FIG. 1A which illustrates a method of enhancing the 
stability of a recursive estimator process. In particular, of a 
process that generates estimates X*(t n + ) of state parameters 
X(t„) of a system from state measurements Y(t„) (wherein t„ 
represents successive times and t n + represents times that 
each lag a respective one of the successive times t n ). 

A first process step 22 of the method generates at least one 
system state measurement Y(t„). With an error covariance 
P(t n ) and a measurement-noise covariance R(t„) of the 
system, process step 23 calculates a gain K(t n ). In process 
step 24, the measurement Y(t„) is then processed with the aid 
of the gain K(t n ) to generate at least one estimate X*(t n + ). 

Process step 25 recursively repeats the generating, calcu- 
lating and processing steps (steps 22-24) to reduce the error 
covariance P(t„) to operational values P 0 that are not sub- 
stantially reduced with further recursion. Finally, process 
step 26 periodically resets the error covariance P(t n ) to a 
reset value P r that exceeds the operational values P 0 to 
thereby reset the calculating step (step 23) and enhance the 
stability of the recursive estimator process. 

An understanding of the processes of FIG. lAis enhanced 
with reference to diagram 110 of FIG. 3 which shows an 
exemplary plot 112 of the error covariance P(t„) as a 
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function of time. The plot 112 illustrates changes imposed 
upon the error covariance P(t n ) by the processes of FIG. 1A. 

In accordance with process step 25 of FIG. 1A, the plot 
112 shows that the error covariance P(t n ) is recursively 
reduced to an operational value 113 that is not substantially 
reduced with further recursion. As a result of the recursive 
estimator process, the error covariance P(t n ) initially falls at 
a rapid rate but that rate declines until the error covariance 
P(t w ) reaches the operational value 113 after which, it is not 
substantially reduced with further recursion that progresses 
along the time axis (e.g., in the time zone 114). 

The operational value 113 represents a steady state value 
for the the error covariance P(t n ) in the absence of state 
perturbations. This steady state, however, is interrupted by 
perturbation inputs that change the state of the system that 
is associated with the processes of FIG. 1A. If the system is 
an attitude control system on a spacecraft, exemplary per- 
turbation inputs are torques applied to the spacecraft’s body 
(e.g., via momentum wheels, thrusters or the movement of 
solar wings). 

The perturbation inputs cause the error covariance P(t n ) to 
temporarily increase to greater perturbation values (as exem- 
plified by the perturbation value 115) before the process step 
25 of FIG. 1A returns it to the steady state operational value 
113. Although the specific magnitudes of the perturbation 
values are functions of the respective system that is associ- 
ated with the processes of FIG. 1A, they will have an an 
upper bound 116. Therefore, the error covariance will have 
an upperound value V ub at the upper bound 116. 

In accordance with process step 26 of FIG. 1A, the error 
covariance P(t n ) is periodically reset to a reset value P r that 
exceeds the operational values. This is exemplified in FIG. 
3 by the reset 118 to a reset value P r which exceeds the 
steady state value 113 and, preferably, exceeds the upper 
bound 116. After the reset 118, the recursive process again 
returns the error covariance to the steady state value 113. 

The relationship between the reset value and the upper- 
bound value V ub may be mathematically expressed as P r = 
(l+a)P Mfc wherein a is a positive number (in applications in 
which the error covariance is represented as a matrix, the 
reset value is formed as the sum of the error covariance 
matrix and a diagonal matrix). Although a is also a function 
of the respective system that is associated with the processes 
of FIG. 1A, simulations indicate that it will generally have 
a value less than 0.1 in spacecraft attitude control applica- 
tions. 

The method 20 of FIG. 1A thus periodically resets the 
error covariance P(t n ) so that the recursive estimator process 
is repetitively forced to start from a selected error covariance 
reset value P r , After each reset, the process continues for a 
time period that is short compared to its total operational 
time period. Because the time period in which the process is 
required to maintain its numerical stability is therefore 
significantly reduced, the demand on its numerical stability 
is also significantly reduced. The method thus not only 
enhances stability but also facilitates the use of less com- 
putationally expensive algorithms. 

The teachings of the invention also facilitate the verifi- 
cation of process stability. This is exemplified by the method 
30 of FIG. IB which is directed to verification of the stability 
of a recursive estimator process that must maintain stability 
for an operational time period T 0 . 

This method begins with the process steps 22-25 of FIG. 
1A. The last of these steps processes the measurement Y(t n ) 
with the aid of the gain K(t n ) to thereby generate at least one 
estimate X*(t n + ) and reduce the error covariance P(t„) to 
operational values P 0 that that are not substantially reduced 


6 

with further recursion. Process step 32 then continues the 
repeating step (process step 25) for at least one reset time 
period T r whose duration is less than that of the operational 
time period T 0 . 

5 At the end of each reset time period T r , the error cova- 
riance P(t n ) is reset in process step 33 to a reset value P r that 
exceeds the operational values P 0 to thereby reset the 
calculating step of process step 24. Finally, step 34 confirms 
the stability of the recursive estimator process over at least 
one reset time period T r . 

In the method 30, the error covariance is reset after each 
reset time period T r . Because the recursive process is 
therefore forced to start from the same covariance after each 
reset time period T r , confirming stability of the process over 
at least one reset time period substantially confirms stability 
15 over the greater operational time period T 0 . 

Because they periodically reset the error covariance to a 
predetermined reset value P r , the methods 20 and 30 of 
FIGS. 1A and IB inhibit numerical instabilities and thus 
maintain a stable recursive estimator process. In a process in 
20 which the error covariance is processed as a matrix, for 
example, the methods insure that the error covariance matrix 
will always be positive definite (a matrix is positive definite 
when for any nonzero vector Z with matching dimensions, 
the product Z r PZ is a positive scalar). 

25 Exemplary applications of the methods of FIGS. 1A and 
IB are in spacecraft attitude control systems such as the 
system 80 of FIG. 4 which controls attitude of a spacecraft 
body 81. In this application, an exemplary operational time 
period T 0 is greater than one year and an exemplary reset 
30 time period T r in the verification process steps 32-34 of FIG. 
IB is on the order of one day. 

The system 80 includes a torque generation system 82, 
gyroscopes 84 and gyroscope data processor 85, star trackers 
86 and star data processor 87, a filter processor 88 and an 
35 attitude control processor 90. The system 82 uses torque - 
generation elements (e.g., thrusters 92, momentum wheels 
93 and/or magnetic torquers 94) to apply torque inputs U(t„) 
to the spacecraft body 82 that change its state X(t„) (e.g., its 
attitude and attitude rate). 

40 The gyroscopes 84 and the star trackers 86 produce 
measurements Y r (t„) and Y(t„) of the state X(t„) from which 
the filter processor 88 generates estimates X*(t n ). Finally, 
the attitude control processor 90 generates attitude control 
signals 92 in response to the attitude difference between the 
45 estimates X*(t„) and a commanded attitude 94. The attitude 
control signals 92 are then applied to the torque generation 
system 82 to generate the torque inputs U(t n ). 

The filter processor 88 is preferably a digital processor 
that is programmed to perform the processes shown in FIG. 
50 2 and the concepts illustrated in the flow charts of FIGS. 1A 
and IB. In an exemplary processor, a star ID (identification) 
and residual computer 100 generates the residue Y(t„)-Y* 
(i n ~) in response to the estimate Y(t„) and inputs from a star 
catalog 102. Again calculator 104 provides the gain K(t w ) (in 
55 accordance with the processes of the gain calculator 44 of 
FIG. 2) which is multiplied in a multiplier 106 to produce the 
correction K(t n ) {Y(t„)-Y*(t„“). 

An attitude and gyro bias propagator 108 uses the rate 
measurement Y r (t„) to extrapolate the estimates X*(t n ) into 
60 the attitude prediction X*(t„“) (that corresponds to a suc- 
cessive measurement) which is summed in a summer 109 
with the correction K(t„) {Y(t„)-Y*(t„ - ) to generate the 
succeeding estimates X*(t n ). In the attitude control system 
80, the resetting step 26 of FIG. 1A would be implemented 
65 in the gain calculator 104. 

FIG. 5 illustrates an exemplary spacecraft 120 that travels 
in an orbital path 122 about the earth 124. Carried in its body 
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81, the spacecraft has a communication system that incudes 
antennas 126 and a power generation system that includes 
solar panels 128. The attitude control system 80 of FIG. 4 is 
preferably included in the spacecraft’s body 81 to control its 
attitude. 5 

For example, the sattelite’s attitude may be controlled to 
direct its antennas 122 at various targets (e.g., a sub -orbital 
point 124 or respective communication cells 126 in a 
terrestrial area 128). The methods of FIGS. 1A and IB are 
implemented in the spacecraft’s processors (e.g., the pro- 
cessor 88 of FIG. 4) to insure the long-term stability of the 
attitude control and to facilitate verification of that stability 
prior to the spacecraft’s insertion into its orbital path 122. 

Although they were described as having matrix form in 
the description of FIG. 2, estimator terms such as the gain 
K(t w ) may, in general, be expressed in any form that facili- 15 
tates their processing. 

The methods of the invention were illustrated with refer- 
ence to Kalman filtering processes but they may be practiced 
with any recursive estimator process. 

The preferred embodiments of the invention described 20 
herein are exemplary and numerous modifications, varia- 
tions and rearrangements can be readily envisioned to 
achieve substantially equivalent results, all of which are 
intended to be embraced within the spirit and scope of the 
invention as defined in the appended claims. 

We claim: 

1. A method of enhancing the stability of a recursive 
estimator process that generates estimates X*(t„ + ) of state 
parameters X(t„) of a system from measurements Y(t n ) 
wherein t n represents successive times and t n + represents 30 
times that each lag a respective one of said successive times 

t n , said method comprising the steps of: 

generating at least one measurement Y(t n ); 

with an error covariance P(t n ) and a measurement-noise 35 
covariance R(t n ) of said process, calculating a gain 
K(t„); 

processing said measurement Y(t n ) with the aid of said 
gain K(t„) to generate at least one estimate X*(t n + ); 

recursively repeating said generating, calculating and 40 
processing steps to reduce said error covariance P(t n ) to 
operational values P 0 that are not substantially reduced 
with further recursion; and 

periodically resetting said error covariance P(t n ) to a reset 
value P r that exceeds said operational values P 0 to 45 
thereby reset said calculating step and enhance said 
stability. 

2. The method of claim 1, wherein said operational values 

P 0 have a steady state value and said reset value P r exceeds 
said steady state value. 50 

3. The method of claim 1, wherein said operational values 
P 0 have an upper bound V ub and said reset value P r exceeds 
said upper bound V ub . 

4. The method of claim 3, wherein said resetting step 
includes the step of determining said reset value P r in 55 
accordance with P r =(l+a)P Mfc wherein a is a positive num- 
ber greater than zero. 

5. The method of claim 1, wherein said at least one 
measurement Y(t n ) comprises a plurality of measurements 
Y(t„) that include a measurement of at least one of said state 60 
parameters X(t n ). 

6. The method of claim 1, wherein said processing step 
includes the steps of: 

with a state transition <D of said process, extrapolating a 
delayed version of an estimate X*(t n _ 1 ) into a state 65 
prediction X*(t n “) and a measurement prediction 


8 

differencing said measurement Y(t„) and said measure- 
ment prediction Y*(t n “) to form a residue Y(t n )-Y* 
(t„-); and 

summing said state prediction X*(t„“) with a correction 
K(t„){Y(tJ-Y*(t n - )} that is the product of said gain 
K(t n ) and said residue Y(t„)-Y*(t„ - ) to thereby gener- 
ate said attitude estimate X*(t n + ); 

wherein t n ~ represents times that each lead a respective 
one of said successive times t„. 

7. The method of claim 6, wherein said extrapolating step 
includes the step of multiplying said state prediction X*(t n “) 
by a measurement converter H(t n ) of said process to realize 
said measurement prediction Y*(t„“). 

8. The method of claim 6, further including the step of 
facilitating said calculating step with said measurement 
converter H(t„) and a process-noise covariance Q(t„) of said 
process. 

9. The method of claim 1, wherein said system is an 
attitude-control system of a spacecraft and said measure- 
ment Y(t n ) includes measurements of at least one of the 
attitude and attitude rate of said spacecraft. 

10. A method of verifying the stability of a recursive 
estimator process for an operational time period T 0 of a 
system wherein said process generates estimates X*(t n + ) of 
state parameters X(t„) of said system from measurements 
Y(t n ) wherein t n represents successive times and t n + repre- 
sents times that each lag a respective one of said successive 
times t„, said method comprising the steps of: 

generating at least one measurement Y(t n ); 

with an error covariance P(t„) and a measurement -noise 
covariance R(t„) of said process, calculating a gain 
K(t„); 

processing said measurement Y(t„) with the aid of said 
gain K(t„) to generate at least one estimate X*(t„ + ); 

recursively repeating said generating, calculating and 
processing steps to reduce said error covariance P(t„) to 
operational values P 0 that are not substantially reduced 
with further recursion; 

continuing said repeating step for at least one reset time 
period T r whose duration is less than that of said 
operational time period T 0 ; 

at the end of each reset time period T r , resetting said error 
covariance P(t n ) to a reset value P r that exceeds said 
operational values P 0 to thereby reset said calculating 
step; and 

confirming stability of said process over at least one reset 
time period T r ; 

said continuing, resetting and confirming steps thereby 
verifying stability over said operational time period T 0 . 

11. The method of claim 10, wherein said operational 
values P 0 have a steady state value and said reset value P r 
exceeds said steady state value. 

12. The method of claim 10, wherein said operational 
values P 0 have an upper bound V ub and said reset value P r 
exceeds said upper bound P ub . 

13. The method of claim 12, wherein said resetting step 
includes the step of determining said reset value P r in 
accordance with P r =(l+a)P Mfc wherein a is a positive num- 
ber greater than zero. 

14. The method of claim 10, wherein said operational time 
period T 0 is greater than one year and said reset time period 
T r is on the order of one day. 

15. The method of claim 10, wherein said system is an 
attitude-control system of a spacecraft and said measure- 
ment Y(t n ) includes measurements of at least one of the 
attitude and attitude rate of said spacecraft. 
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16. The method of claim 10, wherein said processing step 
includes the steps of: 

with a state transition <D of said process, extrapolating a 
delayed version of an estimate X*(t n _ 1 ) into a state 
prediction X*(t n “) and a measurement prediction 5 
Y*(t„-); 

differencing said measurement Y(t n ) and said measure- 
ment prediction Y*(t n “) to form a residue Y(t n )-Y* 
0»'); and lQ 

summing said state prediction X*(t n “) with a correction 
K(t„){Y(t n )-Y*(t n - )} that is the product of said gain 
K(t„) and said residue Y(t n )-Y*(t n “) to thereby gener- 
ate said attitude estimate X*(t n + ); 

wherein t n ~ represents times that each lead a respective 
one of said successive times t n . 

17. The method of claim 16, wherein said extrapolating 

step includes the step of multiplying said state prediction 
X*(t„") by a measurement converter H(t„) of said process to 
realize said measurement prediction Y *0n )* 20 

18. The method of claim 16, further including the step of 
facilitating said calculating step with said measurement 
converter H(t n ) and a process-noise covariance Q(t n ) of said 
process. 

19 . The method of claim 10, wherein said error covariance 2 s 
P(t„) is expressed as a matrix and said confirming step 
includes the step of verifying that said matrix remains a 
positive definite matrix. 

20. The method of claim 10, wherein said system is an 
attitude-control system of a spacecraft and said at least one 30 
measurement Y(t n ) includes measurements of the attitude 
and attitude rate of said spacecraft. 

21. An attitude -controlled spacecraft, comprising: 

a spacecraft body; 

at least one star tracker that is coupled to said body for 35 
providing attitude signals; 

at least one gyroscope that is coupled to said body for 
providing attitude-rate signals; 

a data processor in said body that is programmed to 4Q 
perform a recursive estimator process that generates 
attitude estimates X*(t n + ) from attitude measurements 
Y(t„) wherein t n represents successive times and t n + 
represents times that each lag a respective one of said 
successive times t n , said process including the steps of: 45 

a) generating at least one measurement Y(t„); 

b) with an error covariance P(t n ) and a measurement- 
noise covariance R(t n ) of said process, calculating a 
gain K(t„); 

c) processing said measurement Y(t n ) with the aid of 5Q 
said gain K(t„) to generate at least one estimate 

x*(t„ + ); 

d) recursively repeating said generating, calculating 
and processing steps to reduce said error covariance 
P(t w ) to operational values P 0 that are not substan- 
tially reduced with further recursion; and 
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e) periodically resetting said error covariance P(t„) to a 
reset value P r that exceeds said operational values P 0 
to thereby reset said calculating step; 

an attitude controller in said spacecraft that generates 
torque generation signals in response to an attitude 
difference between a commanded attitude and said 
attitude estimate X*(t n + ); and pi a torque generation 
system that is coupled to generate torques in said body 
in response to said torque generation signals to thereby 
reduce said attitude difference; 

said resetting step thereby resetting said calculating step 
and enhancing said stability. 

22. The spacecraft of claim 21, wherein said operational 
values P 0 have a steady state value and said reset value P r 
exceeds said steady state value. 

23. The method of claim 21, wherein said operational 
values P 0 have an upper bound P ub and said reset value P r 
exceeds said upper bound P ub . 

24. The spacecraft of claim 21, wherein said processing 
step includes the steps of: 

with a state transition O of said process, extrapolating a 
delayed version of an estimate X*(t n _ 1 ) into a state 
prediction X*(t„) and a measurement prediction 

Y*(t„'); 

differencing said measurement Y(t n ) and said measure- 
ment prediction Y*(t„“) to form a residue Y(t„)-Y* 
(t„-); and 

summing said state prediction X*(t„“) with a correction 
K(t„){Y(t„)-Y*(t„ - )} that is the product of said gain 
K(Q and said residue Y(t n )-Y*(t n “) to thereby gener- 
ate said attitude estimate X*(t„ + ); 

wherein t n ~ represents times that each lead a respective 
one of said successive times t„. 

25. The spacecraft of claim 24, wherein said extrapolating 
step includes the step of multiplying said state prediction 
X*(t n ") by a measurement converter H(t„) of said process to 
realize said measurement prediction Y*(t n “). 

26. The spacecraft of claim 24, further including the step 
of facilitating said calculating step with said measurement 
converter H(t„) and a process-noise covariance Q(t„) of said 
process. 

27. The spacecraft of claim 21, wherein said at least one 
star tracker comprises three star trackers directed along three 
mutually -orthogonal axes of said spacecraft and said at least 
one gyroscope comprises three gyroscopes arranged to 
detect spacecraft rotation about said axes. 

28. The spacecraft of claim 21, further including solar 
wings coupled to said body to generate power. 

29. The spacecraft of claim 21, further including antennas 
coupled to said body to transmit and receive communication 
signals. 

30. The spacecraft of claim 21, wherein said torque 
generation system includes thrusters coupled to said body. 





